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ABSTRACT 

We performed J- and 7?-band linear polarimetry with the 4.2 m William Herschel 
Telescope at the Observatorio del Roque de los Muchachos and with the 1.6 m tele- 
scope at the Observatorio do Pico dos Dias, respectively, to derive the magnetic field 
geometry of the diffuse molecular cloud surrounding the embedded protostellar system 
NGC 1333 IRAS 4A. We obtained interstellar polarization data for about two dozen 
stars. The distribution of polarization position angles has low dispersion and suggests 
the existence of an ordered magnetic field component at physical scales larger than the 
protostar. Some of the observed stars present intrinsic polarization and evidence of be- 
ing young stellar objects. The estimated mean orientation of the interstellar magnetic 
field as derived from these data is almost perpendicular to the main direction of the 
magnetic field associated with the dense molecular envelope around IRAS 4A. Since 
the distribution of the CO emission in NGC 1333 indicates that the diffuse molecular 
gas has a multi-layered structure, we suggest that the observed polarization position 
angles are caused by the superposed projection along the line of sight of different 
magnetic field components. 

Subject headings: ISM: clouds - ISM: individual objects: NGC 1333 - ISM: magnetic fields 
- Polarization - Stars: Individual: 2MASS - Techniques: polarimetric 
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1. Introduction 



Infrared and optical polarimetry is a suitable tool for observing magnetic fields within 
molecular clouds at large scales. At these wavelengths p o lariza tion can be produced by dichroic 



extinction of background starlight. 



Davis & Greenstein 



(|l95lh proposed that a fraction of 



non-spherical interstellar dust grains become aligned perpendicular to the local magnetic field due 
to paramagnetic relaxation. Although this mechanisr n is commonly invoked in the literature, it 



seems to be inefficient within molecular clouds (e.g.,|Lazarian 



2007[) . However, a more realistic 



scenario was proposed by several authors who have successfully modeled the perpendicular 
alignment between grains and rnagnet i c fields by radiative torques propelled by anisotropic 



radiation (|Draine & Weingartner 



1996 



Lazarian & Hoang 



20071: 



Hoang & Lazarian 



2008L 



2009 ). 



Aligned dust grains behave like a polarizer to any incoming radiation, absorbing and 
scattering the component of the electric field (E-vectors) parallel to their longest axis. Therefore, 
the observed radiation will carry some degree of linear polarization. The resulting polarization 
map outlines the geometry of the magnetic field lines projected onto the plane of sky (POS). 
Near-infrared (near-IR) polarimetric observations trace visual extinctions of a few tens of 
magnitudes, providing deeper photometry than optical wavelengths. However, the increase in 
interstellar extinction is not usually accompanied by a linear increase in the degree of polarization. 



This has been interpreted as a decrease in the po 



increasing visual extinction ([Goodman et al 



1992 



arization efficiency, or 



19951 



Gerakines et al. 



depo 



1995 



arization, w ith 



Arce et al. 



1998 ) 



Nevertheless, this d epolarization at optical wavelengths is not observed in the Pipe Nebula 



(Franco et al. 



2010h and submillimeter polarization observations show that there is unequivocal 



evidence that grains do align in dense environments with high visual extinction (e.g.. 



Whittet et al 



2008 



Vaillancourt et al.ii2008|) . Additionally, the scattering of stellar light by dust grains also 



generates linear polarization in the optical and near-IR. This type of polarization is found in 
reflection nebulae associated with disks and envelopes of young stars. 
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*^GC 1333 is the most active star-fomiing site in the Perseus molecular cloud (|Lada et al 



1996|) . A large porti on of the NGC 1333 young stellar cluster is composed by low-mass stars 



younger than 1 Myr (|Wilking et al 



20041). In ad dition, there are numer ous embedded pro to stars 



powering molecular and Herbig-Haro outflows ([Knee & Sandell 



20001) ■ There is evidence that 



the rnolecular cloud in NGC 1333 is being distii rbed by the large amount of outflow (jWarin et al 



19961 : 



Sandell & Knee 



2001 



Ouillen et al 



20051). 



and 



The first polarimetric o bservations toward N GC 1333 were carried out by IVrba et al.l (|l976|) 



Tumshek et al, 



( 1980 ) 



Tamura et al, 



(|1988|) conducted ^-band polarimetric observations 



towards the center of the NGC 1333 reflecti on nebula. A larger po larimetric survey covering 



the full Perseus complex was carried out by 



Goodman et al. 



( 19901) . These observations show 



that there is a bimodal distribution of polarization PA., indicating that there are two large scale 
magnetic field components along the line of sight. 

NGC 1333 IRAS 4A (hereafter IRAS 4A), a low-mass protostellar system, has become the 
textbook case of a collapsing magnetized core: high angular submm polarimetric observations 
have revealed that the m agnetic field has an hourglass morphology at scales of few hundred AUs 



( Girart et al. 



1999 



2006b . This is the magnetic field morpholog y predicted by theoretical mode ls 



based on magnetically controlled molecular core collapse (e.g.. iShu et al. 



19871 : 



Mouschovias 



200 ih . Indeed, the synthetic polarization m aps constructed using models of collapsing magnetized 



cores (Galli & Shu 



(|Goncalves et al 



19931 



Shu et al, 



2006h reproduced quite well the observations in IRAS 4A 



20081). In this context, it is worth mentioning that it is still a question of ongoing 



debate whether magnetic fields or in terstellar turbu 



evolution of a molecular cloud (e. g.. 



Crutcher et al. 



€nce plays a major role in the dyn amical 



2009 



Mouschovias & Tassis 



20101) . 



In this paper, we report on one of the first scientific results ob tained with the near-IR ca mera 



LIRIS (Long-slit Inter mediate Resolution Infrared Spectrograph: 



Acosta-Pulido et al. 



2003 



Manchado et al. 



2004|) in its polarimetric mode. The observations were done using the 7-band 
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filter toward stars located relatively close to IRAS 4A (~ 4'-8')- The fields were selected to avoid 
the most active star- forming portion of the NGC 1333 cloud, so the measured polarized light is 
mainly due to dichroic absorption. In order to ascertain the quality of the near-IR data, we also 
provide complementary i?-band linear polarimetry obtained with the Observatorio do Pico dos 
Dias toward the same region. The scientific goal of this work is to compare the magnetic field 
observed in the IRAS 4A molecular core with the larger scale field associated with the cloud 



surrounding IRAS 4A. iGirart et al 



(|2006|) have already done this comparison but using very few 



distant stars (~ 14'-20') retrieved from the iGoodman et al 



(119901) survey. 



2. Observations 

2.1. Near-infrared observations 

The near-infrared observations were carried out in December 2006 and December 2007 at 
the Observatorio del Roque de los Muchachos (La Palma, Canary Islands, Spain). The LIRIS 
camera, attached to the Cassegrain focus of the 4.2 m William Herschel Telescope, is equipped by 
a Hawaii detector of 1024 x 1024 pixels optimized for the 0.8 to 2.5//m range. 

LIRIS is capable of performing polarization observations by using a wedged double 
WoUaston devic e, WeDoWo , which is composed by a combination of two WoUaston prisms and 



two wedges (see 



Oliva 



1997 



1 for detailed description). In this observing mode, the polarized flux 
is measured simultaneously at four different angles (0°, 45°, 90° and 135°). An aperture mask of 
4' X r is used in order to avoid overlapping between the different polarization images. Figured] 
shows a typical LIRIS image in polarimetric mode. The degree of linear polarization can thus be 
determined from data taken at the same time and with the same observing conditions. In order 
to achieve accurate sky subtraction, a 5-point dither pattem was used. Offsets of about 20" were 
adopted along the horizontal, long mask direction. During the 2006 and 2007 campaigns, we took 
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seven and six exposures, respectively, of 20 s per dither position. The 5-point dither cycle was 
repeated several times until completion of the observation. The total observing time for each field 
was 2800 s in 2006 and 2400 s in 2007. 

We carried out 7-band polarization observations of ten fields, six of them with the telescope 
rotator at 0° and four with the rotator at 90° (see Table[I]). Figure |2] indicates the observed fields 
as black and red rectangles, corresponding to observations with rotator at 0° and 90°, respectively. 
We covered the area surveyed by observing with the rotator at 0° and 90°, except for the two upper 
fields. This procedure allows us to compare both data sets and, consequently, to achieve higher 
precision in the estimated polarization parameters. 



2.2. Optical observations 

The optical i?-band linear polarimetry was performed using the 1.6 m telescope of the 
Observatorio do Pico dos Dias (LNA/MCT, Brazil) during observing runs conducted in 2007 
and 2008. A specially adapted CCD camera composed by a half- wave rotating retarder followed 
by a calcite Savart plate and a filter wheel was attached to the focal plane of the telescope. The 
half-wave retarder can be rotated in steps of 22°5 and one polarization modulation cycle is fully 
covered after a complete 90° rotation. The birefringence property of the Savart plate divides 
the incoming light beam into two perpendicularly polarized components: the ordinary and the 
extra-ordinary beams. From the difl'erence in the measured flux for each beam one estimates the 
degree of polarization and its orientation in the plane of the sky. F or a technical descriptio n of 



1 

Magalhaes et all (i 19961) . The 



this polarimetric unit, we refer the interested reader to the work by 
obtained optical data is part of an ongoing large scale (~ 1 square degree) survey whose results 
will be discussed in a forthcoming paper (Franco et al., in preparation). The area covered by 
the optical survey overlaps the portion of the sky observed in near-IR, and in order to make a 
comparative analysis of the results obtained at both wavelengths, we included the optical results 
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gathered for stars lying in the overlapped area in the discussion. 



3. LIRIS data reduction and calibration 



The near-IR data reduction was performed using the lirisdr package developed by the LIRIS 
team in the IRAF environmentQ Given the particular geometry of the frames (see Fig.[D, the 
first procedure was to slice the image into four frames. Each set of frames corresponding to a 
given polarization stage is processed independently. The data reduction process comprises sky 
subtraction, flat-fielding, geometrical distortion correction, and finally co-addition of images after 
registering. A second background subtraction was performed upon flat-fielded images in order to 
avoid the residuals introduc ed by the vertical gradient due to the reset anomaly eff"ect associated 



with the Hawaii arrays (e.g.. 



Acosta-Pulido et al. 



2006|) . An approximate astrometric solution was 



determined based on the image header parameters. 



3.1. Photometry 

Aperture photometry of the field stars in each slice was obtained using the task Object 
Detection, available within Starlink Gaia software^ The aperture radius used was ~ 4", which 
corresponds to 3 times the median seeing of the night. The background was extracted from an 
annulus with an inner radius of 6" and an outer radius of 8". The astrometric solution of each 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by 
the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

^GAIA is a derivative of the Skycat catalogue and image display tool, developed as a part of 
the VLT project at ESO. 
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slice was tweaked using the astrometric tools available within the Starlink Gaia software. We used 
the 2MASS catalogue to perfomi the photometric and astrometric calibrations. In our sample, we 
reached J magnitudes as faint as ~ 17. As a final step, we identified the counterparts of each object 
in the four slices in order to compute the polarization properties. In some cases, matching of stars 
observed with rotator at 0° and 90° was also necessary since some objects were present in both 
sets of observations. 



3.2. Polarimetric analysis 

Using the WeDoWo, we measured simultaneously four polarization states in each of the 
strips as, 

1 

(2) 
(3) 
(4) 

where /*, 2* and are the Stokes parameters of the object to be measured, and the factors 
^[0,90,45, 135] represent the transmission for each polarization state. In this case, the normalized 
Stokes parameters can be determined by. 



io - 190 ^0/90 

q* = ■ ^ ■ , — (5) 

Iq + Z90 ro/90 

''45 ~ '135 ^45/135 

u* = : — — ; (6) 

'45 + '135 M5/135 

where the factors fo/9o and ^45/135 measure the relative transmission of the ordinary and 
extraordinary rays for each WoUaston. These factors were calibrated using non-polarized 



ioiPA 


= 0] = 


= -?o(I* + Q*) 




= 0] = 


= 2 (I* ~ Q*) 




= 0] = 


1 

= 2^45(I*+U,) 


h35[PA 


= 0] = 


1 

- 2 ^135 (I* ~ U*) 
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standards and resulted in the values fo/9o = 0.997 and ^45/135 = 1.030, with an uncertainty of about 
0.002 in both cases. 

The rotation of the whole instrument by 90° causes the exchange of the optical paths for the 
orthogonal polarization vectors. Now, the resulting polarization states are given by 



ioiPA 


= 90] = 






= 90] = 


1 




= 90] = 




il35[PA 


= 90] = 





(7) 
(8) 
(9) 
(10) 



This effect can be used in order to get a more accurate estimate of the Stokes parameters 
because the combination of both measurements, PA=0° and 90°, results in the cancelation of the 
transmission factors and reduces flat-field uncertainties. The normalized Stokes parameters are 
then computed by 



= R^V '^"g = kiPA = 90]/i,o[PA = 90] ^''^ 

= g^i^,being4= -•--[^^ = Q]/-"-[^^ = QJ (12) 
Ru + l ^ ^ h5[PA = 90]/iu5[PA = 90] ^ ' 

Finally, after estimation of the q and u Stokes parameter, the degree of linear polarization and 
the position of polarization angle (measured eastwards with respect to the North Celestial Pole) 
are calculated as 



= ^^ql + ul (13) 



1 _i /m 



e = -tanM-|. 



Flux errors in iq, zgo, Z45 and Z135 are dominated by photon shot noise while the theoretical 
error in polarization fraction was estimated performing error propagation through the previous 
equations. In addition, we calculated the errors in p using a Monte Carlo method, which returned 
values similar to those estimated from error propagation . The Icr uncertainty in 6 was estimated 



(z) by applying the relation derived by 



Serkowski 



1974|) using standard error propagation, that 



is, CTfl = 28°. 65 (Tp/p, when p / cTp > 5; or (z'z) graphically with the aid of the curve proposed by 
Naghizadeh-Khouei & Clarkd (119931) when p/cfp < 5. 



Figure |3] shows the polarization uncertainty as a function of the i-band magnitude achieved 
with our LIRIS observations. The observed distribution suggests that the uncertainties are 
dominated by photon shot noise, as expected for a sample collected with fixed exposure time. As 
expected, the uncertainties decrease when the data taken at 0° and 90° are combined. There is 
a natural limit which is due to the uncertainty bias when measuring low levels of polarization. 
Bias in the degree of linear polarization (p) comes from the fact that this quantity is defined 
as a quadratic sum of q and u, which produces a non-zero polarization esti mate due to the 



uncer t ainties in their measureme nt (for a detailed discussion see for instance. ISimmons & Stewart 



1985 



Wardle & Kronberg 



1974|) . In order to remove th e polarization bias a n d com pute the 



true polarization, we used the prescription proposed by ISimmons & Stewart! (|1985|) for low 
polarization stars. The true polarization degree can be approximated by the expressions ptme = 
if PobJcTp < Ka, otherw ise p,r„p = ipi},. - crj ■ KlYI ^. We adopted Ka = I, which corresponds to 



the estimator defined by 



Wardle & Kronberg 



( 1974 ). 



3.3. Standard stars 



Observations of polarized and unpolarized standard stars were taken in order to calibrate 
the instrumental characteristics of LIRIS in its polarimetric mode. Table |2] summarizes the 
general information for these stars: columns 1 to 8 indicate their name, equatorial coordinates. 
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type, polarization degree and position angle, the filter used for the polarization measurements 
and the reference, respectively. Unpolarized standard stars were observed to check for any 
possible instrumental polarization and for systematic errors in our polarimetry. The unpolarized 
stars G191B2B and BD+28d4211 were observed with rotator at 0° and 90°. The polarized 
intensity measured for the two unpolarized standards was very small (see Table |3]): The measured 
normalized Stokes q and u were 0.051% and 0.226%, respectively, with the rotator at 0°, 
and -0.117% and 0.119%, respectively, with the rotator at 90°. Table [3] shows the observed 
polarization degree before and after bias correction for the two unpolarized standards. The 
measurements taken in the two epochs for BD+2 8d421 1 give consistent values. We applied the 



method proposed by 



Simmons & Stewart (1985) for a 99% confidence level for the observed 



unpolarized standards. This resulted in a small, if any, instrumental polarization. 

The polarized standard star CMa Rl No. 24 was observed in order to verify the zero 
point of the polarization position angles. Table |4] summarizes the results obtained for the four 
measurements conducted for this object. As expected, high quality data are less sensitive to 
biasing, and the unbiased polarization has basically the same values of the observed polarization. 
Ta king into account the uncertainties, we see that our y-band data match the result obtained 



by 



Whittet et al. 



(Il992l) . The differ ence between the av erage P.A. obtained for these four 



measurements and that obtained by 



Whittet et al. 



(Il992h is ~ 6°, which is very close to the 



statistical deviation of our measurements thus discarding any further correction for the zero angle 
calibration. 
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4. Polarization properties 
4.1. Infrared data 

Table |5] contains a summary of our near-IR and optical polarization data for stars with 
a signal-to-noise in the polarization intensity higher than unity. Column 1 gives the star's 
identification number in our catalogue. Columns 2 and 3 show the equatorial coordinates. Column 
4 gives the 7-band magnitude. Columns 5 to 12 show the polarization degree and the polarization 
P. A. (with their uncertainties) for the R and J bands. The last two columns indicate the rotator 
position used to acquire the near-IR data, and the object type, respectively. Figure |4] shows that, 
excluding star 13, the polarization RA. distribution measured in the 7-band is quite narrow with 
a mean position angle of 160° and a standard deviation of only 12°. We note that the 7-band 
polarization uncertainties may be overestimated. First, the aforementioned standard deviation is 
about half the mean l-cr uncertainty of the polarization RA., and second, there is a very good 
agreement between the near-IR and optical data (see next section). 

Figure|5] shows the spatial distribution of the near-IR polarization vectors overlaid on 
the 2MASS 7-band image. The polarized stars with a declination below 31° 10' have larger 
polarization degrees than those above this value. This subsample comprises most of stars with 
mean RA. ^ 160°. Star number 13 is the only object in our catalogue presenting intrinsic 
polarization (see section |5]). 



4.2. Comparison with optical data 



Previous optical polarimetric observations performed toward the field of vi ew shown in Fig 



detected only two polari zed stars of the 7-band sample, stars number 2 and 13 (IVrba et al 



Menard & Bastien 



1976 



19921), and our data are in good agreement with them. The 7?-band polarimetric 



sample has 12 stars in common with our near-IR data. Figure |6] shows the polarization vectors 
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in both bands plotted over a DSS image. There is noticeably good agreement between the two 
polarization data sets (see also Fig. |7]). Thus, the mean value of the P.A. difference between the R 
and 7-band for the 12 stars is 6°5. 

Figure[8]represents the resulting Fnir vs. Fvisitie diagram for the 12 stars with both J and R 
polarimetric measurements. The wavelength at which the polarization is highest, A,„ax, is related 
to the mean size of the interstellar grain s responsible for producing the observed polarized light 



( Serkowski et al, 



1975 



McMillan 



interstellar medium is 0.55 fim (jSerkowski et al. 



19781). The typical value of A,nax observed for the diffuse 



1975h. However, in molecular clouds A, 



1990 



Whittet et al. 



2007). In the case of 



appears to change with the visual extinction (|Andersson & Potter 
it has been shown that A„r,x changes between 0.66 and 0.89 /im (iTumshek et al. 



NGC 1333, 



19801: 



Cemis 



1992!). The two solid lines in Fig.[8]show the expected relation for A,„ax = 0-55 



and 0.86 fim. The data do not show any preferential regime of maximum polarization due to the 
low statistics. Polarization observations over a wider range of wavelengths are necessary to refine 
this characterization. 



5. Polarization from YSOs and foreground stars 



Our near-IR and optical maps are characterized by a uniform component predominant to the 
south of the IRAS 4A/4B double system (see Fig. [6l). However, north of IRAS 4A/4B the few 
detected polarized stars have a broader angle distribution. Most of these stars are likely young 
stellar objects (YSOs). In these cases, the polarization is produced by intrinsic scattering within 



circumstellar disks rather than by interstellar abso rption. Several authors have studied the physica l 

2 



Brown & McLean 



1977 



Mundt & Friec 



1983) 



properties of YSOs by means of polarimetry (e.g., 
and, in general, observations suggest that near-IR polarization vectors, when produced by single 
scattering, are oriented perpendicularly to optically thin disks, while multiple scattering within 



optically thick disks generates polarization vectors whose P.A. are parallel to their long axis 



-14- 



(lAngelll 19691 : iBastien & Menardlll990l : IPerevra et al.ll2009h 



The YSOs possibly showing intrinsic polarization in the near-IR and/or 7?-band are listed 
below: 



LkHcf 271 (star n. 13 of Table|5]l is a Classical T Tauri star (jLada et al.lll974|) . The near-IR 
polarization angle and degree are in excellent agreement with the i?-band data. Previous 
observations showed that the polarization varies considerably, which has being interpreted 
as arising from outburst s or inhomogeneities in a circumste liar shell within an optically 



thick circumstellar disk (Tamuraetal 



1988 



Menard & Bastien 



1992h 



• SVS 13A (star n. 24) was observed only in the i?-band. The obtained polarization is in 
excellent agree ment with the value p reviously measured in the ^-band (Py = 7.2 + 0.9% an d 



P. A. = 56 + 4°, 



2002 



Tamura et al. 



Anglada et al 



2004 



1988). 



Chen et al. 



his ob ject is a well-studied source (IRodriguez et al. 



3) 



20091) which powers a bi polar and collimated outflo w 



assoc iated with the well-known Herbig-Haro objects HH 7-11 (jHerbig 



1974 



Strom et al 



1974|) . The orientation of the outflow is roughly perpendicular to the polarization PA., 



which suggests that the disk is optically thick. 

• 2MASS J03290289+31 16010 (star n. 23) is close to SVS 13A and has a very low degree of 
polarization. According to the SIMBAD Astronomical Database, this is a bright (J ^ 12.8) 
K-type star and it is possibly a foreground star. In fact, previous spectral an alysis and 



photometric studie s place this star at a distance of only ~50 pc from the Sun (|Aspin et al. 



1994 



Aspin 



20031) . 



• ASR 8 (star n. 25) is classified as a brown dwarf by SIMBAD. However, an extensive 
survey on the evolutio nary state of s tars in NGC 1333 identifies this object as a T Tauri star 



with a mass of 0.7 M p, (lAspin 



(IGetman et al. 



20031) . which is reinforced by the presence of X-ray emission 



20021). We therefore attribute the optical polarization measured for this star 
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due to intrinsic scattering. 

In addition, there are two bright infrared stars (J < 13.0) with low polarization (stars n. 28 
and 34 in Table[5]) that are apparently not associated with YSOs, as no star formation or nebulosity 
signs has been reported in the literature. Their 2MASS color indices suggest that they may be 
unreddened M-type dwarf stars, which is also corroborated by the low degree of polarization. We 
therefore consider these objects to be foreground stars. 



6. The magnetic field in NGC 1333 



6.1. The distribution of dust and molecular gas in NGC 1333 



The most detailed picture of the distri bution of gas and dust in the Perseus cloud has 



been provided by the COMPLETE project (|Ridge et al 



2006al : 



Pineda et al, 



20081), a survey of 



near/far-infrared extinction data, and of atomic, molecular, and thermal dust continuum emission 
obtained over a large area. These data show a wide range of visual magnitudes for NGC 1333, 
and a non-Gaussian CO spectral profile consistent with multi-velocity components. These results 
are consiste nt and likely related to a layered cloud structure along the line of sight, which was first 



proposed by 



Ungerechts & Thaddeus 



( 1987 ). Interstellar extinction studies of field stars toward 



NGC 1333 also suggest a t least two components in the line of sight at diff'erent distances toward 



NGC 1333 (Cemis 



1990 ) 



2006bh . the region 



According to column density maps of the Perseus cloud (|Ridge et al.l 
studied here lies in the lower density envelope of NGC 1333. Maps of high density molecular 
tracers (NiH^, HCO^) as well as of the 870 fim dust emission, show that around IRAS 4A the 
dense gas has a fila mentary distribution oriented in the N\y-SE direction, wit h the long axis 



positioned at -142° (|Sandell & Knee 



2001 



Olmi et al 



2005 



Walsh et al 



2007h 
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6.2. The field morphology as traced by the diffuse gas 



The near-IR and optical polarization vectors of the background stars shown in Fig. |6] trace 
the POS component of the magnetic field associated with the lower density envelope around 
IRAS 4A/4B. South of these sources, where we have most of the polarization sample, the 
magnetic field has a direction of - 160°. The observed confi guration is consisten t with the results 



obtained at much larger s cale by 



Goodman et al, 



the COMPLETE survey dRidge et al 



( 1990 ) and 



Tamura et al, 



(Il988h . According to 



2006ah the polarization was measured toward regions with a 



visual extinction of 4 to 5 mag. 



The magnetic field orientation derived fro m our data is 



filamentary structure associated with IRAS 4A (IWalsh et al. 



rough 



2007 



y parallel to the dens e 



Sandell & Knee 



20011) . 



However, the submm polarization maps towards IRAS 4A and IRAS 4B show that t he magnetic 



field \ yithin the filament is app roximately perpendicular to the filament's major axis (|Girart et al. 



1999 . 



2006 



Attard et al 



20091), and is therefore perpendicular to the magnetic field direction 



traced by our optical and near-IR data. The single-dish submm polarization map from 

2) 



Attard et al. 



(120091) around IRAS 4A is associated with visual extinctions as low as ~10 magnitudes, which is 
a typical value for near-IR extinction data. Therefore, the submm and near-IR/optical data seem 
to reveal substantial changes in the magnetic field topology between the dense filament and the 
diffuse molecular envelope that surrounds it. Such a sharp twist in the field is hard to explain 
by means of structural changes in the magnetic field only, because within the observed field, the 
position angle of the optical and near-IR polarimetric data is quite uniform (see Fig. ID. Instead, 
the two data sets may be simply tracing distinct gas components. As explained in § 16. 1[ there 
is observational evidence of a multi-component structure for the NGC 1333 molecular cloud. 
Figure |5] shows the ^^CO and ^^CO spectra extracted from a box containing the region studied 
here. These spectra show at least three distinguishable velocity components: a faint emission 
centered at vlsr - 2 km s^' (seen more clearly in the '^CO data), the peak of the '^CO data 
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centered at ~7.6 km s and the peak o f the CO d ata at ~6.7 km s . This last component has the 



same vlsr of the IRAS 4 A dense core (|Choi 



200 ih . Therefore, whereas the submm polarization 



measurements trace only the molecular cloud component associated with the IRAS 4A dense 
core, the near-IR and optical polarimetric data are probably tracing the mean magnetic field of the 
different velocity molecular cloud components observed in the CO maps. Nevertheless, further 
observations are needed in order to obtain a more complete description of the magnetic field in 
this region. 



7. Conclusions 

We have carried out one of the first polarimetric observations in the 7-band collected 
with the WHT/LIRIS infrared camera. We also present 7? -band linear polarimetry obtained at 
the Observatorio do Pico dos Dias. We observed an area of ~ 6' x 4' around the NGC 1333 
IRAS 4A/4B protostellar system. The main conclusions of this work are: 

• The infrared polarization map derived for the surveyed area is highly consistent with the 
optical map obtained with a different telescope and observational technique. Therefore, 
the near-IR polarimetric capabilities of LIRIS have proved to be scientifically trustworthy 
for the astronomical community, and assure this mode will be useful for gathering 
measurements of objects experiencing high interstellar extinction inaccessible to optical 
instruments. 

• The polarization map obtained for the surveyed area is dominated by a well-ordered 
component produced by dichroic interstellar absorption. However, there are objects, some 
of them catalogued as YSOs, that show a transversal component which may be generated 
by intemal scattering within circumstellar disks. 

• The magnetic field morphology traced by the near-IR/optical map is almost perpendicular 
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with respect to the field morphology obtained with the submillimeter data toward the dense 
molecular core around IRAS 4A/4B. The near-IR/optical polarimetric data trace the field 
morphology of the diffuse molecular gas, which is known to have a multi-velocity structure. 
That is, the observed resulting magnetic field direction is probably the averaged magnetic 
field over several distinct velocity components of the cloud. CO molecular data obtained for 
this line of sight show non-Gaussian line profiles that are consistent with this hypothesis. 
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Table 1 . Log of the observations^ 



Target 


0^2000 






Obs. date 


Rotator 


ID 


(hh:mm:ss.ss) 


(dd:mm:ss.ss) 




n 


Fl 


03:29:22.01 


+31 


:09:42.12 


2006 Dec 26 





F2 


03:29:25.24 


+31 


:08:41.88 


2006 Dec 26 





F3 


03:29:23.98 


+31 


:07:49.41 


2006 Dec 26 





F4 


03:29:22.90 


+31 


:06:54.62 


2006 Dec 26 





F5 


03:29:34.26 


+31 


: 12:50.08 


2007 Dec 13 





F6 


03:29:20.52 


+31 


:15:59.63 


2007 Dec 13 





Flp 


03:29:15.44 


+31 


:08:12.82 


2006 Dec 26 


90 


F2p 


03:29:19.59 


+31 


:08:28.20 


2007 Dec 13 


90 


F3p 


03:29:25.70 


+31 


:08:17.24 


2007 Dec 13 


90 


F4p 


03:29:30.95 


+31 


:08:30.44 


2007 Dec 13 


90 



''The night of 2006 December 27 had very limited weather con- 
ditions and only calibrators were observed. 
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Table 2. Standard stars 



ID 


Q'2000 


52000 


Type 


P 


6" 


Filter 


Ref. 




(hh:mm:ss.sss) 


(dd:mm:ss.ss) 




(%) 


n 






CMa Rl No. 24'' 


07:04:47.364 


-10:56:17.44 


Polaiized 


2.1 ±0.05 


86 ± 1 


J 


1 


BD+28d4211 


21:51:11.070 


+28:51:51.80 


Unpolarized 


0.041 ± 0.031 


38.66 


N*^ 


2 










0.067 ± 0.023 


135.00 


U 


2 










0.063 ± 0.023 


30.30 


B 


2 










0.054 ± 0.027 


54.22 


V 


2 


G191B2B 


05:05:30.621 


+52:49:51.97 


Unpolarized 


0.065 ± 0.038 


91.75 


u 


2 










0.090 ± 0.048 


156.82 


B 


2 










0.061 ± 0.038 


147.65 


V 


2 



"Position angles measured from north to east. 

''This is the only polaiized standard stai' with well established polarization propeities in the /-band. 

'^ "Near-UV" filter cen tered in 3450 A and with full width at half maximum (FWHM) bandpass of 650 A. For details, 
see lSchmidt et alj j 19921) . 



References. — (l) IWhittet et alj )1992|) : (2) ISchmidt et al.l )1992|) , the two unpolarized stars are taken from a com- 
pilation of optical calibration data collected with the Hubble Space Telescope. 



Table 3. Observational results for the unpolarized standard stars 



ID 


Mission 


^obs 


Via 


P 

^ true 


P /P " 

^ mml ^ max 






(%) 






(%) 


G191B2B 


2007 


0.41 


1.38 


0.28 


0.00/1.10 


BD+28d4211 


2006 


0.09 


0.60 


0.00 


0.00/0.41 


BD+28d4211 


2007 


0.15 


1.17 


0.08 


0.00/0.45 



''Minimum and maximum va lue for the degree of polar - 
ization at 99% confidence level (ISimmons & Stew ar] |l985i . 



Table 4. Observational results for the polarized standard star 



(%) (%) (%) (%) C) C) 



CMaRl 2.10 0.26 7.98 2.08 1.38/2.78 92.7 4 

No 24 2.04 0.26 7.99 2.02 1.35/2.70 94.3 4 

2.33 0.26 9.06 2.31 1.62/3.00 93.5 3 

2.05 0.22 9.43 2.04 1.44/2.62 86.8 3 



^Minimum and maxim um value for the degree of p olarization 
at 99% confidence level ( Simmons & Stewartlll985 ). 
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TableS. Polarization data 



ID 


^^2000 


"^2000 


7" 


V 






P 




ftme 




(Tu 
^ f/ 


Rotator position 


Class'' 






^ LIVJ.. JiiJii . A A. f 


(mag) 








1'%') 


1'%') 




\ ) 


V / 






1 


03:29:14.58 


31:06:38.20 


14.70 


4.04 


0.75 


151.5 


1.74 


1.27 


1.41 


173 


27 


0,90 


IR source 


2 


03:29:14.89 


31:09:27.50 


12.67 


4.53 


0.06 


157.8 


3.10 


0.64 


3.03 


157 


7 


0,90 


Star 


3 


03:29:15.05 


31:08:06.90 


16.06 








3.16 


2.43 


2.02 


160 


33 


0,90 


IR source 


4 


03:29:16.08 


31:07:31.40 


12.99 


4.51 


0.12 


166.2 


2.32 


0.74 


2.20 


172 


10 


0,90 


IR source 


5 


03:29:16.20 


31:07:34.00 


13.83 


3.92 


0.21 


157.6 


1.88 


0.97 


1.61 


167 


18 


0,90 




6 


03:29:16.68 


31:16:18.30 


10.74 


0.87 


0.14 


117.5 


1.67 


0.60 


1.56 


135 


11 





IR source 


7 


03:29:17.52 


31:07:33.20 


15.04 


3.88 


0.25 


163.3 


1.82 


1.45 


1.10 


171 


37 


0,90 


IR source 


8 


03:29:17.91 


31:07:07.70 


15.34 








2.70 


1.91 


1.90 


169 


29 


0,90 


IR source 


9 


03:29:18.21 


31:07:55.70 


14.50 


3.32 


0.48 


167.2 


2.55 


1.43 


2.11 


163 


20 


0,90 


IR source 


10 


03:29:18.64 


31:09:59.60 


12.50 


4.40 


0.04 


156.4 


4.68 


0.65 


4.63 


155 


4 


0,90 


Star 


11 


03:29:20.01 


31:09:54.30 


12.45 


5.75 


0.10 


164.0 


4.89 


0.64 


4.85 


160 


4 


0,90 


Star 


12 


03:29:20.10 


31:08:54.00 


16.08 








3.15 


2.44 


1.99 


155 


33 


0,90 




13 


03:29:21.87 


31:15:36.30 


11.33 


0.94 


0.09 


49.4 


1.60 


0.75 


1.41 


49 


16 





CTTS'= 


14 


03:29:23.50 


31:07:25.00 


16.57 








4.54 


3.46 


2.94 


167 


33 


0,90 




15 


03:29:24.70 


31.07:27.00 


17.16 








4.56 


4.42 


1.12 


167 


24 


0,90 




16 


03:29:25.60 


31:08:43.00 


17.00 








4.90 


3.99 


2.85 


172 


23 


0,90 


IR source 


17 


03:29:27.04 


31:08:04.60 


12.39 


4.61 


0.29 


157.6 


3.53 


0.66 


3.47 


169 


6 


0,90 


IR source 


18 


03:29:27.16 


31:06:48.20 


14.92 


5.20 


0.51 


166.1 


3.32 


1.79 


2.80 


173 


20 


0,90 


IR source 


19 


03:29:28.99 


31:10:00.30 


13.36 








4.01 


0.93 


3.90 


142 


7 


0,90 


IR source 


20 


03:29:29.60 


31:08:47.90 


16.30 








3.17 


2.71 


1.64 


156 


42 


0,90 


IR source 


21 


03:29:30.80 


31.06:33.00 


17.97 








7.01 


6.51 


2.61 


160 


25 


0,90 


IR source 


22 


03:29:32.41 


31:13:01.10 


13.34 








1.90 


1.21 


1.47 


135 


24 





IR source 


stars with i?-band data only 


23 


03:29:02.87 


31:16:00.82 


12.84 


0.68 


0.15 


66.4 














YSOC 


24 


03:29:03.74 


31:16:03.60 


11.74 


7.58 


0.54 


57.7 














V5I2Per 


25 


03:29:04.04 


31:17:06.66 


13.31 


1.38 


0.52 


86.3 














BDf 


26 


03:29:07.39 


31:10:49.02 


13.10 


3.11 


0.17 


167.0 














Star 


27 


03:29:09.57 


31:09:08.68 


14.94 


4.61 


0.61 


161.1 














Star 


28 


03:29:12.16 


31:08:10.91 


12.99 


0.39 


0.08 


71.8 














IR source 


29 


03:29:17.84 


31:05:37.40 


14.04 


4.35 


0.45 


164.4 














IR source 
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Table 5 — Continued 



ID 


O'2000 


(52000 






o-pR 


(fib 




Ptrue 




(Tg^ Rotator position 


Class'' 




(hh:mm:ss.ss) 


(dd:mm:ss.ss) 


(mag) 


(%) 


(%) 


n 


(%) (%) 


(%) 


n 


D (°) 




30 


03:29:20.59 


31:06:11.48 


15.28 


5.13 


0.93 


174.5 










IR source 


31 


03:29:29.11 


31:06:08.77 


14.07 


5.36 


0.32 


167.3 










IR source 


32 


03:29:34.24 


31:07:53.33 


13.64 


2.15 


0.39 


153.7 










IR source 


33 


03:29:39.77 


31:14:51.61 




1.38 


0.19 


13.6 












34 


03:29:40.43 


31:12:46.38 


13.04 


0.34 


0.06 


41.4 










IR source 



^IRIS J magnitude for stars observed with this camera, 2MASS / magnitude for stars observed only in the R-band. 
''Position angles are counted from north to east. 

'^1(T uncertainty of the position angle (see text for explanation on how it was estimated), 
''object's type as found in the SIMBAD Astronomical Database. 



'=Classical T-Tauri Star l lLada et al.ll974 . 



'Brown Dwarf jitoh et aLlboiol) 




Fig. 1. — A typical CCD image in polarization mode. The four strips correspond to the 0°, 90°, 
135° and 45° polarization vectors from which the Stokes parameters are calculated. 
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Fig. 2. — DSS 7? -band image of our science targets. Black boxes are observed fields with rotator 
at 0° while red boxes are observed fields at 90°. Crosses mark the positions of the protostars NGC 
1333 IRAS 4A and NGC 1333 IRAS 4B. The ellipse in the upper right comer indicates the star- 
forming region, where no science targets were select in order to avoid polarization data due to dust 
scattering. 
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Fig. 3. — Distribution of the polarization uncertainty (crp) with respect to the J magnitude obtained 
for each field star with the telescope rotator at 0° (green triangles), 90° (red circles) and a combi- 
nation of both (crosses). This plot contains only stars whose signal-to-noise ratio of the combined 
setup (crosses) is better than 1 (except for three stars that also have a signal-to-noise ratio better 
than 1 but which were only observed with the rotator at 0°, see Table |5]). Note that uncertainties 
are lower when the combination of images taken with the telescope rotator at 0° and 90° is used. 
The large discrepancy observed between 0° and 90° errors for some stars are due to the distinct 
observation epochs of each data set. 
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Fig. 4. — Distribution of polarization angles of the near-IR data. The histogram is binned in 10° 
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Fig. 5. — y-band polarization vectors in NGC 1333 plotted over a 2MASS y-band image. Vector 
length scale is shown on the upper left comer. Green vectors indicate stars with P/crp > 3 while 
red vectors have 1 < P/crp < 3. Open circles indicate positions of observed objects with P/crp < 1. 
Some of the detected polarized stars in the /-band does not have a 2MASS counterpart, indicating 
that the obtained LIRIS data probe deeper visual extinctions in the cloud than the 2MASS. 



-34- 



5% 



15 



o 
o 
o 

c 
o 

ca 
_c 

"(5 
Q 



31° 10' 



Qgh 45S 



X IRAS 4 A 
TrAS 4B 




30'' 1 5' 

Right Ascension (2000) 



00^ 



28"^ 45^ 



Fig. 6. — Comparison between optical (blue vectors) and near-IR (red vectors) data. The polari- 
metric map is plotted over a DSS 7?-band image. The vector length scale is shown on the upper 
left comer. Orange vectors represent the averaged magnetic field of IRAS 4A and IRAS 4B, as 



obtained by submillimeter observations of 



Attard et al. 



(20091). 
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Fig. 7. — Comparative diagram of the position angles obtained for the optical (i?-band) and near-IR 
(7-band) data sets. 
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Fig. 8. — Spectral Energy Distribution of the observed linear polarization in near-IR and visible. 
Solid lines indicate constant A,„ax of 0.55 and 0.86 ftm from bottom to top, respectively. Dashed 
lines represent constant p^ax of 2, 4 and 6% from the origin going outwards, respectively. 
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Fig. 9. — Averaged spectra of the ^^CO 1-0 (solid line) and the ^''CO 1-0 (dashed line) lines 
obtained over a region of about 5' centered at the position a(J2000)=3''29'"24'^and 5(J2000)=3 1° 8'. 



This region covers the F1-F4 and Flp-F4p observed fiel ds with the WHT (see § 12.11). The CO 



spectra was retrieved from the COMPLETE data archive (jRidge et al. 



2006al: 



The dotted vertical line shows the systemic velocity of the IRAS 4A core (jChoi 



Pineda et al.1 



20081). 



2001b. 



